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Abstract: 

Extensions of free relativistic particle models, living in Non-Commutative space-times that might be 
relevant in the context of Quantum Gravity in Planck scale physics are well known. In the present 
work the dynamics of such charged particles in the presence of external electromagnetic interactions are 
studied in detail. Furthermore we have also considered the electromagnetic field to be dynamical and 
have derived the modified forms of Lienard-Wiechert like potentials for these extended charged particle 
models. In all the above cases we exploit the new form of K-Minkowski algebra where electromagnetic 
effects are incorporated in the lowest order, in the Dirac framework of Hamiltonian constraint analysis. 



I. INTRODUCTION 

From a theoretical perspective, it appears that the existing theoretical framework of physical laws might 
require a significant modification in the regime of Quantum Gravity that is for Planck scale physics. A 
particularly interesting scenario, termed as Doubly Special Relativity (DSR), was pioneered by Amelino- 
Camelia [1] . (For reviews see for example [2]). This model is an extension of Einstein's Special Relativity. 
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In Special Relativity there is one (inertial) observer independent dimensional parameter, the velocity of 
light but in DSR there are two observer independent dimensional parameters, the velocity of light and 
an energy scale, conveniently chosen as the Plank energy. Since it is expected that Planck scale provides 
the borderline of qualitatively different physics, (related to the nature of spacetime in particular), it is 
intuitively obvious that all inertial observers should agree on the same value of the Planck scale. The 
mathematical framework of DSR appears most natural in Non-Commutative (NC) geometry [2, 3] and 
a well studied NC phase space is the one obeying ^-Minkowski NC algebra [4]. The NC geometry was 
introduced in physics very early by [5] and resurrected to its present interest by the work of [6]. Out 
of several alternative bases of the algebra [3], that are inter-related through non- linear transformations, 
we will concentrate on the Magueijo-Smolin (MS) base [7]. In general, extended relativistic principles in 
DSR requires a modification in the Einstein energy-momentum relation for a free particle of rest mass mo 
and momentum P M , P^P^ = m,Q. Quite naturally the canonical Lagrangian L — m^yj {(dX ^) / (dr)) 2 also 
needs to be drastically altered to account for the DSR particle dispersion relation and /{-Minkowski NC 
phase space. Two of the present authors [8], along with many other workers [9], have provided explicit 
models for DSR particles that is consistent with the /{-Minkowski NC algebra and DSR dispersion relation. 
Many novel features of the DSR particle have come in to light from the analysis of these models. 

It needs to be emphasized that mostly these works deal with the free DSR particle but many interesting 
properties as well as drawbacks of the model (if there are any) are not encountered unless the particle 
is subjected to non-trivial interactions. In the present paper we have precisely done that. We have 
considered a DSR particle interacting with a U(l) gauge theory. We have mostly concentrated on treating 
the gauge field as external but we have also commented on the DSR effects on the Lienard-Wichert type 
potentials arising from the charged DSR particle. 

The paper is organized as follows: in Section II we introduce the /{-Minkowski NC algebra in MS base 
and the interacting DSR particle - U(l) gauge theory. In Section III we derive the extended /{-Minkowski 
algebra in presence of U(l) gauge interaction and obtain the equations of motion. Section IV deals with 
explicit solutions for particle trajectories in a perturbative framework that provide valuable insights of 
the DSR effects on a charged particle motion. In Section V we comment on the electromagnetic potentials 
induced by the charged DSR particle. The paper concludes in Section VI with discussions and future 
outlook. 

II. INTERACTING CHARGED MS PARTICLE IN U{1) FIELD 

The /{-Minkowski NC phase space algebra is defined as, 

{x\x a } = - ■ {x\x*} = ; {x\p>} = -gV ; {p*p"} = 0. (I) 
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Our metric is diag g 00 = —g lt — 1 and k is the NC parameter. Rest of the phase space algebra is given 
below, 



{x°,p i }=p i /K; {x\p°} = 0; {x°, P } = -l+p°/K. 
The above is rewritten in a covariant form, 



(2) 



(3) 



where rj = l,rji = 0. This algebra appeared in [10]. Detailed studies of similar types of algebra are 
provided in [3]. This algebra has emerged before in an earlier work of one of us [11] where it was 
embedded in a more general algebra. For n —¥ oo one recovers the normal canonical phase space. Using 
K-Minkowski NC algebra (3), it is straightforward to check that the Lorentz generators — x^p^ — x^p^ 
satisfy the canonical Lorentz algebra J^ v , 

{J^ 7 J a/3 } = g^J va + g^jP" + g ul3 J a v + g va . (4) 

However the particle dispersion relation that remains invariant under (3) with respect to J^ iV treated as 
rotation generators in the NC framework is given below: 



P 



1 



(VP) 



1 



E 

K 



(5) 



We use the shorthand notation (ab) = a^b^. This is the well known MS dispersion relation that reduces 
to the Special Theory relation for k — > oo. 

In order to derive a Lagrangian for the DSR particle ab initio it is quite difficult to guess how to 
proceed since it has to be compatible with the complicated K-Minkowski symplectic structure (3) . In this 
perspective, it is very convenient to exploit the Darboux-like canonical variables. The idea is to construct 
the Lagrangian in terms of canonical variables which are in fact the Darboux map of the physical non- 
canonical (or NC) degrees of freedom. The next step is to reexpress the the Lagrangian in term of the 
physical NC variables via the inverse Darboux map. We have successfully exploited this technique in a 
previous work [8]. In the present case, with the canonical phase space, 



{X^tPv} — —g^u] {X^jX^} — {P^,P V } — 0, 



the inverse Darboux map is 



X„ 



(vp) 



(jml) (i-f) 



(6) 



(7) 
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whereas the Darboux map is, 

,..X.(l + ^).X.(l + ^) t y „ = 7 -L = II |^. (8) 



K 



We consider the first order form of a particle interacting minimally with U(l) gauge field, 
L = P"{t)X^t) - \(P 2 (t) - ml) + ( d'ye S(y - X(r))X^(y) 



P»(t)X,(t) - A(P 2 (r) - m 2 ) + e X,(t)A»(X(t)) 



(px)(r)p) X 2 2 f (r)p)\ \ A UIV \( ■ fi ivp)\ (VP) 



In the last equality we have used the inverse Darboux map (7) and the notation a^b^ = (ab). As has 
been shown before [8], for the free theory e = 0, Dirac constraint analysis [12] yields the NC K-Minkowski 
symplectic structure. 

III. NONCOMMUTATIVE SPACE AND DYNAMICS OF INTERACTING MS PARTICLE 

IN EXTERNAL [7(1) FIELD 

Our first task is to compute the symplectic structure modified by the electromagnetic interaction. We 
utilize the Dirac formalism [12] for Hamiltonian analysis of constraint system. A similar analysis for the 
free theory, (e = 0) was done in [8]. From (II) the constraints are, 

^ = n--p M -e^«o, (10) 



<-i? + ^ + ^«o, (11) 

where ip = 1 — and the conjugate momenta are 

IP- — -v + ethA ■ IP- dL - {PX)V » < xA ^ hoi 
The constraints are non-commuting since 

(Ik -f£) <13 » 
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in 2 n 2 \ — ^BHH X v r l^) 



(15) 



Hence the constraints are Second Class in Dirac sense and we wish to compute the Dirac bracket [12], as 
defined below, 



{A, B}* = {A, B} - {A, rfHrf, if}~ 1 {7f, B}, 
so that the constraints can be strongly put to zero. Expressing the constraint matrix as 





and using the inverse 



a + 



■% F uXX X V» 



= A + eB 



(A + eB)- 1 = A' 1 - e(A~ 1 BA~ 1 ) + 0(e 2 ) 
the required inverse of the constraint matrix to O(e) is computed below: 



9u\ - ^mPu - ^x\r] a F av ei/j 2 [F u \ + ^{p^Fax ~ p\F av )\ 



Therefore to 0(e) the Dirac Brackets are 



{x a ,Xp}* = {x a ,xp} - {x a ,ril}{r]l,r]{} 1 {r){,Xp} = -(x a rjp - XpT} a ) 

K 



(16) 



(17) 



(18) 



(19) 



1 eijj 2 

{x a ,PfiY = {x a ,P{s} - {x a ,r] v }{n v ,T) x }~ {VxiPp} = (-9ap + -VaPp) H x a r} v F v p (20) 

{p a ,P[sY = {Pa,Pp} - {Pa,r) v }{ r )v,'nx}~ {Vx>Pp} = # ( F <*p H (p a F v p - ppF va )). (21) 

For e = the above reduces to the free ^-Minkowski algebra [3]. It is worth pointing out that {x ai xp} 
remains unchanged and the modifications involve the gauge invariant field tensor F fll/ (X). 

To obtain the equations of motion we write the Hamiltonian first in terms of canonical variables and 
subsequently in terms of NC variables: 



(22) 



As P 2 — Toq => p 2 = rrtQi/j 2 . Note that this H should not be identified as the energy but it is easy to 
check that in the canonical sector it clearly reproduces the correct Hamiltonian equations of motion, 



(dx a )/dT) = x a = {x a ,HY, p a = {p a ,H}* 



(23) 
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We should mention that there is an approximation involved in deriving H in (22). Strictly speaking we 
should have used the Darboux map pertaining to the U(l) modified algebra (19 - 21) whereas, mainly for 
simplicity we have used the Darboux map (7,8) for the free NC algebra (3) as derived in [8]. Although 
0(e) correction to the Darboux map is straightforward to compute we have not considered this here. 
But as we demonstrate in rest of the paper, even in this relatively simplified situation, the K-corrections 
generate strikingly qualitative differences. 

Therefore the equations of motion to 0(e) follow: 



x a = (2 - V) 



1 

mo 



1 



(i]Fp)x c 



P 



(2-VO; 



m 



F al3 pp + -{r}Fp)p a + 2mlF afi r ] p 



(24) 



(25) 



Let us write down the modified Lorentz force equation to order 0(e). We differentiate (24) with respect 
to t to obtain, 



x a = (2 - V) 



1 

too 



1 



(r)Fp)x c 



2ip 



moip 2 



(26) 



Now using (24) we replace all p-dependent terms in the R.H.S. of (26) and keep only 0(e) terms throught 
This excercise yields, 



m x a = eF a/3 p p + 



(4 - yoy> 4 TOo 

(2-V) 



(rjFx) 



ip 2 mo 



(r/x) -ljx a - tJj m Q (7jFx)x a 
(27) ' 

Clearly the MS Lorentz force law is quite involved with non-linear terms coming in through the k- 
correction and for k — > oo the conventional Lorentz law is smoothly recovered. 



IV. TRAJECTORIES OF CHARGED MS PARTICLES IN EXTERNAL (7(1) FIELD 

We will try to solve perturbatively the first order r-derivative force equation (25) for some specific 

external field configuration. Due to the nature of the basic (MS form of K-Minkowski algebra in (3)) NC 

structure considered, it is easy to see that the K-corrections always involve the electric field in the form 

(rjF) a , as is clear from (19-21). Hence for situations with external purely magnetic fields the dynamics 

does not change very much in a qualitative way. Hence for the present case, we consider only external 

constant electric field. Again this analysis is true for combined electric and magnetic fields but with the 

condition | E \ > | B \ > such that one can Lorentz transform to another frame with electric field only (see 

for example [13]). Without loss of generality, let us consider the electric field to be along the direction of 

xi-axis, Ei = E, E% = E$ = 0. Breaking up the set of equations in (25), we find 

dp eE 2eE 

~T = — Pi H PoPi , 

dT TO, TJIK 
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dpi 
dr 



eE eE , o „ o 9x 
—P0 + (p 2 + 2m 2 -pj), 



dp2 
dr 



eE dp 3 

P1P2 , -r- 

mn dr 



eE 



-PxPs- 



(28) 



We immediately notice a qualitative difference: For the normal particle the force is impressed only along 
the direction of the external electric field and the energy changes accordingly. As for example in the 
present case the electric field is along X\ and for the normal particle the momentum will change only x\ 
direction. However, for the MS particle we observe that not only are the rates of changes energy and X\- 
momentum modified, but more interestingly momenta along 22 and X3 directions are also affected. Note 
that the mass m will actually depend on the initial velocity as we specify below. We further differentiate 
the above set of equation and keeping terms up to the first order term of ^, a little rearrangement leads 
to: 

d 2 po 



dr 2 



1 



Pa + - (3pg + p\ + 2m 2 ) 



(29) 



d 2 pi 
dr 2 



Pi 



-P0P1 



(30) 



<fp2 
dr 2 



-P0P2 



(31) 



d 2 P3 



-P0P3 



(32) 



dr 2 k 
where — = e. 

This is the system of non-linear ordinary differential equation we plan to solve to the leading order in 
1/k. Let the initial conditions for the momenta at r = be po(0) — m, P2(0) — fnvo, Pi(0) — ^3(0) = 
so that m = 7== , v = 7=^7^ where Vq — ^\t=o, vq — ^:\ T =o- We provide some details of solving 



V^ 2 



7^ 



perturbatively the non-linear set of equations (29-32) in the Appendix. To 0(1/ k) the solutions are, 

^ , < , 2 f™ 2 2\ ^ , / s 2 ™ 2 „ wo ^ m 2 + 2m 2 

Po = m Cosh (er) H — — — h m Cosh (er) H — - — Cosh (2erJ , 

K V 6 / 3k k 



2 ( in 2 \ nr 2 

Pi = m Sinh (er) H — ( h m 2 I Sink (er) H Sink (2er), 

k V 6 / 3k 



p 2 = mv 



1 + — (1 - Cosh{eT)) 



, P3 = 0. 



(33) 
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Since p = and pi = m^p-, integrating p^ from the above set of equation with respect to r gives 



, to 2 2m 2 \ , , . m 2 , (to 2 + 2m 2 )er 

tticxq = m+ — 1 binnier) + ——binhlZeT) , 

3k /c / 3k k 



(34) 



to 2 2m 2 \ , to 2 / to 2to 2 

tticxi = m + — 1 Cosh(eT) + -—uosMZct) — to + — 1 

3k k J ok \ 2k K 



(35) 



meX2 = mvo 



er I H binnier) 

k J K 



(36) 



x 3 =0. 



(37) 



Hence we have found the paths of the MS particle in proper time. Finally we can eliminate r from (35) 
and (36) to construct the trajectories of the particle as, 



TO / 5TO\ 1 2 In , 

* 1 + M 1+ 2^j} "if 11 



I \ in \ x\ e 2 E 2 to 2 / 'mi 
Vq 3tok e 2 ^ 2 



(38) 



2k ) v,\ :>!iii; e*£J* \ k 

The /t-corrected hyperbolic path is recovered provided we neglect the x\ term as it is of 0{e 2 /k). This 
leads to a modified hyperbola: 

2 e 2 E 2 



eE ( 5m 
—xx + 1 + — 
m \ Zk 



1 



IIto 
2k 



1 



5m 



The eccentricity S is reduced by the K-correction, 



9 Utow 2 
£ = \/l + ug- 



2 k 



(39) 



(40) 



Various features of the paths and trajectories of the MS particle are graphically shown in Section VI. 



V. LIENARD-WIECHERT POTENTIAL FOR MS PARTICLE 

The full interacting action for electrodynamics of the conventional particle is 



A = J [-(dr)^(A>(r)X M (r)] - J d% e 5(y - X(r))X^{y) -\J*y F^(y)F^(y) 

= J [-(dr)^(X^r)X^T)} - J (dr) e X,(r)A^X(r)) - \ J ' d*y F^(y)F^(y). (41) 

Now we notice that the point particle and interaction sector depend on which are the dynamical 
variables. Indeed, when we consider the MS particle, X^ need to be replaced by the physical z M variables 
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that obey the NC algebra. On the other hand, the argument y M in the pure Maxwell sector are just 
parameters and not dynamical variables which obviously are given by A^(y). Hence in the free Maxwell 
theory and resulting Maxwell equations there is no change due to the MS particle. The only change that 
will appear is in the interacting system via the source term in the Maxwell equation, 

d*Fp, = j„(x), j v = ev v (42) 

Now we can consider the Lienert-Wiechert potential for a moving charged MS particle. For the normal 
free charged particle this retarded potential is obtained from the inhomogeneous equation 

d v d v A» = f, (43) 

in Lorentz gauge d^A^ — 0. This leads to the general solution 

A»(x) = A%(x) + [ dx'G(x - x')f{x'), (44) 



where Aq(x) satisfies the homogeneous equation 9 2 Ag(a;) = (see for example [14]). Using the retarded 
form of the Green's Function G(x — x') one finds the potential for a current , 

A>(x) = A»{x) + IJ ^"^ij/p^'W )- (45) 
The form of current used by us is 

f{x) = e[d i y5{y~x(T))v^y). (46) 



Hence, recalling the argument presented above, for the retarded potential for the MS particle we need 
only introduce the form of MS current which in turn means we should use the form of velocity appropriate 
for the MS particle. 

We return to our equations of motion valid to 0(1/ k): 

dx a n . „ (VP) s P a , e / r? s a 
— 1 + 3 h (riFp)x , 

CLT K TTlQ TOqK 



cLt mo k moK 

Again, to 0(e) and 0(1/k), the above equations can be inverted to read 



(47) 



^-[(1-3^)^ + !^]. (48 ) 
m k cLt k or 

This motivates us to replace the velocity (dx a )/(d,T) in the charge current by the MS velocity 

^^[(l-sMj^ + V^]. (49) 
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First of all we consider the case with no external field, = so that the replacement is simply 

dx a , W^dx a 

17 ^-s^' (50) 

where the energy W is computed from the MS dispersion relation (5), 

W = + J(<2 + (1 _< ( ^ +TO 2 )]/(1 _ Hiy (51) 

K V K K Z K 

To 0{k) this simplifies to 

W~JjP + m 2 Q ) - ^ (52) 



with the rest energy being 



Wb~(l-^)mo. 



Hence, for a MS particle at rest, the factor (1 — 3-^-) ~ (1 — 3—) modifies the Coulomb potential, 

A = - e -(l-3 n ^); A t = 0. (53) 

Defining R = x — j7(to), n = R/R,v = dx/dt, we find the situation is much more involved for a moving 
MS particle because now the factor (1 — 3^) is no longer a constant but depends on velocity, 

e W -> ev W 

4itR(1 — n.v) K 4itR(1 — n.v) k 

However, probably more studies are needed to interpret our results for the MS particle induced retarded 
potential in the presence of external field F^ v due to the presence of the x Q -dependent term in (49) in 
the expression for velocity. We will not pursue this analysis any further in the present publication. 



VI. DISCUSSION, CONCLUSION AND FUTURE OUTLOOK 



We will mainly study the results obtained in Section IV regarding the paths and trajectories of MS 
particles in external electromagnetic field. Let us demonstrate pictorially the K-effects by way of diagrams. 
For the sake of comparison we have always provided a diagram for 1/k = depicting the conventional 
particle. In the Figure 1 we plot the particle energy as a function of r. To our approximation of 0(1/ k) 
the diagrams are not qualitatively different from the 1/r = (conventional particle) case but the energy 
increases more rapidly as 1/k increases. Figure 2 shows the plot of x\ vs. r that is the coordinate in 
the direction of the external field. Once again due to the approximation involved all the curves show 
hyperbolic nature with the graphs with larger 1/k more steeply curved. However in Figure 3, where we 
plot X2 vs. r, there appear qualitatively different paths for the MS particle in relation to the conventional 
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particle. The latter is a simple straight line as there is no acceleration along X2 in the conventional particle 
case but for MS particle the hyperbolic nature dominates for later r due to the K-dependent acceleration 
term in x%- 

In the second set of figures (Figures 4,5 and 6) we show the trajectories in the x\ — X2 plane for 
different choices of external parameters. In Figure 4, we plot the trajectories for different values of k 
for fixed external conditions: constant values of e ~ electric field E in xi-direction and constant vq ~ 
initial energy of the particle. In Figure 5 we keep vq unchanged and consider the effect of varying e on 
different K-valued particles. Figure 6 shows the effect of different w on different K-valued particles where 
throughout e is unaltered. Some common features of the trajectories are the following: Quite expectedly 
the effects are more pronounced as 1/k increases, as is apparent from Figure 4. In both figures 5 and 
6, we have compared the behavior of 1/k = 0.1 particle with 1/k = 0, the conventional particle. We 
notice that the /{-effects are more appreciable for larger values of the external parameters. This is also 
understandable since the external parameters are directly connected to the energy scales of the system 
and all ^-corrections involve the ratio of (energy)/ k. Lastly a general feature is that in all instances of 
fixed e and vq the K-particle trajectory bent more toward the a;2-axis with respect to the conventional 
K = particle. This is because ^-corrections produce a non-trivial acceleration in the ^-direction even 
though the electric field is in xi-direction, in contrast to the conventional particle for which the velocity 
in xi direction remains constant. From the relation(40) also it is clear that the eccentricity £ of the 
hyperbolic trajectory decreases as 1/k increases. 

Finally let us summarize our work. We have studied the behavior of an extended model of charged 
particle, pertaining to the Magueijo-Smolin form of particle in the Doubly Special Relativity framework. 
The free particle resides in a phase space with K-Minkowski form of Non-Commutative symplectic struc- 
ture. As we have argued in the paper, studies of the above free particle model carried out so far in 
the literature is indeed interesting and important but the free particle can not by itself provide a true 
picture. It is imperative to study the particle in the presence of non-trivial interactions as the interacting 
model can show the novelties as well as possible drawbacks (if there are any) of the new exotic particle 
models. In the present paper we have studied in detail the charged Magueijo-Smolin particle interacting 
with electromagnetic field, both external as well as dynamical, in a perturbative framework. As we have 
demonstrated there are non-trivial K-corrections in the particle trajectories. We have briefly discussed 
the dynamical gauge fields as well and how the Lienard-Wiechert type of potentials get modified due to 
the K-corrected source term. Deeper analysis is required in the latter area. 

For future work we intend to take up dynamics of the gauge field more thoroughly in the K-Minkowski 
framework, especially in 2 + 1-dimcnsions that is an area of recent interest [15]. It is interesting to note 
that 2 + 1-dimensional spacetime allows more general forms gauge kinetio term other than just the 
Maxwell term, such as the Chern-Simons term or a combination of both the Maxwell and Chern-Simons 



term (that is the topologically massive gauge theory [16]). We will report on the MS particle dyn 
in these forms of gauge kinetic terms in future publications. 






FIG. 3: Plot of X2 vs. r for fixed e — 0.1, m = 1 and different k with - = for broken line. 
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FIG. 4: Trajectory plot for different k with e = 0.001 ? initial velocity vo = 0.001. 
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FIG. 5: Trajectory plot for fixed vq = 0.01 and different e with ~ = 0.1 (unbroken line) and ^ = (broken line). 
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FIG. 6: Trajectory plot for fixed e = 0.001 and different vq with ^ = 0.2 (unbroken line) and ^ = (broken line). 
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VII. APPENDIX 



We provide some computational details for the solutions of the system of non- linear ordinary differential 
equation (29-32), 



(fpo 
dr 2 



1 



Po + - {3p 2 + pj + 2m 2 ) 



<Ppi 

' d T 2 



Pi + -PoPi 

K 



(55) 



d 2 P 2 
dT 2 



.2 ,ri 
~PaP2, 



< - d 2 p 3 e 2 



dr 2 



-P0P3 



(56) 



We solve it in a perturbative framework in the first non-trivial order in 1/k. Let us assume 
Poo, Pio, P20, P30 to be the solution of po, pi, p2, P3, for k — > oo, that is the conventional case. 
Now we expand p^ about the conventional solutions p^o '■ 



c d g h 

PO = POO + -POI, Pi = PlO + -Pll, P2 = P20 + -P21, P-i = P30 + -P31 

he Hi hh hh 



(57) 



Here c, d, g, h are constants to be evaluated later. If at t = the initial velocity v is along y-axis then 
initial condition for the above set of perturbation are 

p (0) = m =*> poo(0) = m , poi(O) = , 
pi(0) =0 => pio(0) =pn(0) =0 , 
p 2 (0) = mv =j> p 2 o(0) = mv , p 2 i(0) = , 



p 3 (0) = p 3Q (0) =p 3 i(0) = 

Using (57) the above set of differential equation becomes 

d 2 Poo _ 2 d 2 Pio 2 d 2 p 2 o _ „ d 2 p 30 



dr 2 



e Poo, 



dT 2 



£ Pio, 



dT 2 



0. 



dT 2 



d 2 poi _ 2 
Z dr 2 ~" 



[cpoi + 3poo + p\o + 2m2 ] > 



d 2 P2i 2 , d 2 P3i 

- = -e P00P20, h——- = -e 



dT 2 



dT 2 



P00P30- 



(58) 

(59) 
(60) 
(61) 
(62) 
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The solution of (59) are respectively 

Poo = aie eT + a 2 e~ eT , p w = bie eT + b 2 e~ er 



(63) 



P20 = ci + c 2 t, p 30 = d-L + d 2 r 



(64) 



Now using (63) and (64) the solution of the differential equation (60), (61) and (62) are respectively 

3c 3c c 



. r , 2a l & l 2er , 2a 2 fe 2 _ 2er 2(ai& 2 + 6ia 2 ) 

Pn = mie + m 2 e H — ^-^e H e 



3D 



3d 



(66) 



P21 = 



9 L 



2aic 2 . 
aici h raic 2 ) e 



2a 2 c 2 

a 2 ci H h ra 2 c 2 J < 



(67) 



P31 



t - ~ j aid 2 e tT + (t + - ) (i-ulu 



in i mi 
C T + C 



(68) 



Therefore 



po=[ai+ ck\ ^ + ( + ch\ e _ er + 3a? + 6i e2eT + 3a|+5| e _ 2er _ 6^ + 26^ + 2^ 
' K J \ K I 3k 3k ck 



Pi = Ul + 



2ai6i_ e2£r + 2a 2 6 2 c _ 2er _ 2{a 1 b 2 + b 1 a 2 ) 



3 k 



3 k 



P2 = ci + c 2 r 



aici + raic 2 I e 



2a 2 c 2 " 

a 2 Ci + ra 2 c 2 H I e 



c gr eg 



K K 



(71) 



P3 = - 



1 



r - ^ j aid 2 e er + ( r + ? ) </.//,( 



K K 



(72) 



All the arbitrary constants introduced above are solved by checking consistency of the set of solutions 
leading to the result given below: 

m m m 2 2 m 2 2 

ai = a 2 = — , bi = -o 2 = — , cli — cl 2 = — + m , ami = -dm 2 — — + m , 
2 2 o 6 



Cl 77Z 

ci = touq , c 2 = , c' = , c" = -y- , di = , d 2 = , c'" = , c"" = (73) 
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